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We derive the generalized nonlinear wave equations in frequency domain for 
^ two coupled pulses with different polarizations interacting in quadratic nonlin- 

c/^ ear crystals. The equations include for the first time both quadratic and cubic 

• ^ nonlinearities, delayed Raman effects and anisotropic nonlinearities. The full 

E>^ potential of these wave equations is demonstrated by investigating simulations 

p I of sohtons generated in the process of ultrafast cascaded second-harmonic 

generation. We show that a balance in the sohton delay can be achieved due 
^ to competition between self-steepening, Raman effects and self-steep ening-like 

^ effects from cascading originating in the group-velocity mismatch between 

^ the pump and second harmonic. Wc analyze the first-order contributions, 



and show that this balance can be broken to create fast or slow pulses. 
Through further simulations we demonstrate few-cycle compressed solitons in 
extremely short crystals, where spectral phenomena such as blue/red shifting, 
t> non-local radiation and dispersive-wave generation are observed and marked, 

^ which help improving the experimental knowledge of cascading nonlinear 

soliton pulse compression. @ 2012 Optical Society of America 

OCIS codes: 190.2620, 190.5650, 190.7110. 



1. Introduction 

Cascading nonlinearities are known to be produced accompanied with harmonic generation 
in the limit of large phase-mismatch. As an example phase-mismatched second-harmonic 
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generation (SHG) gives rise to a well-known cascaded quadratic nonlinearity which is Kerr- 
like and can be written as a nonlinear refractive index riciL i.e. it induces a nonlinear index 
change that is proportional to the pump intensity. While the energy conversion to the second- 
harmonic is weak, a Kerr-like nonlinear phase shift is induced on the fundamental wave (FW), 
which can be large and its sign can be tuned by the phase-mismatch; importantly a negative- 
sign self-defocusing nonlinearity is accessible. In this case, soliton pulse compression can be 
fulfilled by combining the self-defocusing cascading nonlinearity with normal dispersion: this 
so-called cascading nonlinear soliton pulse compressor operates within the visible and near- 
infrared region and without a power limit compared with the self-focusing pulse compression. 
Hence, few-cycle and high energy soliton pulses can be generated [2]-[4]. 

Nonlinear crystals with low dispersion, such as beta-barium borate (BBO), potassium 
titanyl phosphate (KTP) and lithium niobate (LN), are good candidates for cascading non- 
linear soliton pulse compression as they have a decent 2°*^ order nonlinear index, and large 
cascading nonlinearities can be achieved through tuning the crystal close to phase-matching 
by either exploiting the birefringent (Type-I) interaction between the FW and SH [s] , or by 
exploiting quasi-phase-matching (QPM) technology. The goal is to overcome the material 
self-focusing Kerr nonlinearity of the crystal by the self-defocusing cascading nonlinearity, 
and it can usually be done if the phase-mismatch is taken low enough. Noncritical (Type-0) 
interactions, where FW and SH are polarized along the crystal axes is desirable because it 
exploits the largest 2"^^ order tensor component, but since the process is not phase-matchable 
at all, QPM has historically been used to reduce the residual phase-mismatch and thereby 
increase the cascading strength [6 -10 . In general these nonlinear crystals also promise a 
stable and compact pulse compressor as they have a short soliton length and a large energy 
threshold. 

Theoretical analysis and experimental demonstrations have been widely reported on cas- 
cading nonlinear soliton pulse compression. An elegant theory is that, in the limit of large 
phase-mismatch, harmonics are all considered as perturbations on the FW. Therefore, the 
multiple coupled-wave equations (CWEs) can be reduced to a single one governing the FW 
which is similar to the nonlinear Schroodinger equation (NLSE) [iT 13 . In the NLS-like equa- 
tion of SHG process, the cascading quadratic nonlinearity is scaled by both the quadratic non- 
linearity squared and the inverse of the phase-mismatch. Its properties are revealed through 
a non-instantaneous ("nonlocal") cascading response which has a Lorentzian shape in fre- 
quency domain. This nonlocal response reveals that there are two different characteristic 
regimes for the response function, known as stationary and non-stationary regimes and the 
transition between the two regimes is determined largely by the phase-mismatch parameter 
and the group-velocity mismatch (GVM) parameter 14 . A detailed explanation was made by 
Bache et al. to help understanding most experiments 15 , 16 . For example, the few-cycle high 
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energy pulse compression in Type-I BBO crystal at 1250 nm reported by Moses and Wise [s], 
was conducted in a regime in which the phase-mismatch configuration was in the stationary 
region so that the cascaded response is broadband 14 and allow for few-cycle compressed 
pulses to form. On the other hand, Ashihara et al. reported the moderate pulse compression 
also using Type-1 BBO crystal but launching a pump at 800 nm [2], which is actually in 
the non-stationary region 15 and the input pulses contained a low soliton order. In both of 
their experiments, slow pulses with shock fronts are observed which is induced by GVM as 
1*^* order cascading response. In a latter report by Ashihara et al., Type-I periodically poled 
magnesium-oxide-doped lithium niobate (PPMgOLN) was used to compress the pulse at 
1560 nm where the GVM between the FW and the SH is suppressed |8|. Therefore, there is 
no GVM-related cascading response and the compressed pulse is symmetrical. Recently, the 
first Type-0 few-cycle pulse compression experiment was demonstrated by Zhou et al., which 
used a Type-0 LN bulk crystal pumped at 1300 nm |4], in which the phase- mismatch is lo- 
cated well into the stationary region, see Fig. [l]^b). The compression was accomplished within 
an extremely short length of 1 mm as the pulses contained a high soliton order caused by the 
high intensity and because the FW group velocity dispersion (GVD) in LN is very large. As 
LN has a dominant delayed Raman response (caused by a strong IR phonon coupling) which 
continuously transfers energy from high-frequency (blue) components to low-frequency (red) 
components, tendencies of early pulse splitting were observed. 

In this paper, we investigate theoretically and numerically the self-steepening performance 
and fast /slow pulse tendencies in cascading nonlinear soliton pulse compression induced by 
both the nonlocal cascading response and the material Raman effects. We do that by directly 
solving the forward Maxwell equation (FME), and a counterbalance of fast/slow pulses is 
demonstrated. First, we extend the FME to describe interaction in a birefringent material, 
and model the electric field of the optical wave by directly solving a group of nonlinear 
wave equations in frequency domain (NWEFs), in which both the 2^^^^ order and 3''*^ order 
nonlinear-induced polarizations are included, see Sec|2j From the NWEFs, we present to our 
knowledge the most complete NLS-like equation which includes the system intrinsic self- 
steepening, a full edition nonlocal cascading response and the material Raman effects, see 
Sec|3| Then, we give a review on the conditions of cascading nonlinear soliton pulse compres- 
sion (Secjij) and based on them we discuss the consequences of employing QPM technology. 
Self-steepening performance and the counterbalance of fast/slow pulses are highlighted as 
they are greatly dependent on the phase mismatch parameter, see SecjSj In Secj6| we give 
realistic and convincing simulations to show the self-steepening performance and fast/slow 
pulses dominated by both the nonlocal cascading response and the material Raman effects. 
Since the NWEFs are generalized, typical phenomena such as the self phase modulation 
(SPM) induced broadening in spectrum; generation of dispersive wave and resonant SH 
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radiation in non-stationary regime are still observed. These phenomena can be marked in 
spectrum, which helps improving the spectral knowledge of cascading nonlinear soliton pulse 
compression. Finally, we present our conclusions in SecjT] 

2. Nonlinear Wave Equations in Frequency Domain (NWEFs) 

In this part, we will present our master equation which directly deals with the electric field. 
Starting from the 1+lD wave equation of the electric field and making the slowly evolving 



wave approximation (SEWA), we obtain the FME 17,18 : 



The "~" mark indicates that both the electric field E and the nonlinear induced 
polarization Pnl are written in frequency domain (we use the Fourier transform 
-^[/(^)] ~ I-^ f{t)e~'^^*dt and the forward definition as J^E ^ — 2k(c(;)E). In a uniaxial crys- 
tal, ordinary and extraordinary waves, Eo and E^., are two separate components of the electric 

~ 1/2 

field vector, and the propagation constant, k(ci;) = (w/c) (1 + , is a function of the lin- 

ear electric susceptibility xi.^) > which is a diagonal matrix containing only the self-response 
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Among all types of nonlinear induced polarizations, 2'^'^ order and 3'''^ order nonlinear 
polarization are the most relevant. Normally, the 2"*^ order nonlinear polarization, which 
appears only in non-centrosymmetric media, is considered instantaneous and, in frequency 
domain, its component can be written as: 

where j,ai,a2 = o,e, xf]xia2 effective 2'^'^ order nonlinear susceptibility (with unit 

m/V), which is actually the element of a 3-rank tensor. The 3'''^ order nonlinear polarization 
consists of both the electronic Kerr response (which is modeled as an instantaneous response) 
and the delayed vibrational Raman response. Its component can be written as [20]: 



i(3) 
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(3) 

Xi;Qia203 effective 3'^'^ order nonlinear susceptibility (with unit m^/V^) coming from a 

4-rank tensor, /r indicates the relative fraction of the material Raman effects, h-[i{Lj) is the 
normalized Raman frequency response function, typically modeled as a complex Lorentzian. 
Note that if there is no Raman fraction in the material (which is the case for a noble gas), 
the 3'^'^ order nonlinear polarization will be completely instantaneous. 
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Thus, the FME can be expanded to two equations corresponding to the ordinary and 
extraordinary wave, named nonhnear wave equations in frequency domain (NWEFs): 



dz 



+ iko(uj)Eo = -i 
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(5) 



We note that in order to solve the above equations, the contents of propagation constants in 
negative frequencies should be known. According to the causality and reality of the material 
response, kj{u) 
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The contents in 



to) shows a property of conjugate symmetry 
negative frequencies are therefore linked to that in positive frequencies. Besides, parameters 
such as k{u}), x^'^\ x^^^ ^-^e dependent on the rotation angle {9, Lp) of the crystal reference 



frame 21,22 26 



Obviously, Eq.(|4]) and Eq.([5]) describe the dynamics of the electric field rather than the 
field's envelope. They can be solved by using split step Fourier method together with the 
Runge-Kutta method. The NWEFs, which for the first time include polarization coupling, 
second and third order anisotropic nonlinearities and delayed Raman nonlinear effects are 
quite general in concept Since any of the possible phase-matching conditions (Type-0, -I and 
-II) are easily modeled, and by including the anisotropy of the nonlinearities, the different 
wave-mixing possibilities (3- and 4-wave mixing, including third-harmonic generation and 
parametric up- and down-conversion) are automatically included properly. 

3. Full edition NLS-like Equation 

From our NWEFs, an NLS-like equation governing the FW envelope can also be derived. 
First, we degrade the NWEFs to the common CWEs by employing more approximations such 
as slowly- varying envelope approximation (SVEA), assumption of constant refractive index 
around the carrier frequency and assumption of no overlap between the harmonics. Then, 
in the strong cascading limit A/cL ^ 1 (A/c = k2 — 2ki is the phase-mismatch, subscripts 
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1 and 2 correspond to the FW and the SH, respectively, and L is the interaction length) 
and in the assumption of undepleted FW (i.e. only the nonlinear phase shift is accumulated 
as cascading nonlinearity) , the ansatz U2{z,t) = A2{t) exp{—iAkz) is used and the cross 
phase modulation (XPM) is neglected. Thus, in the FW reference frame we get the NLS-like 
equation which is the most complete compared with what Ilday et al. [l2], Menyuk et al. fl3], 
Moses et al. pT|, Valiulis et al. 23 and Bache et al. 14,15 got because the equation includes 
a full edition nonlocal cascading response as well as the material Raman effects and intrinsic 
self-steepening. The dimensionless form 16 of the complete NLS-like equation is: 



4 - a) U,^.,n(Ak)N!^ (l - Ulir) (l - jy^M 

POO 

(1 - fn) |f/i(r)|'f/i(r) + /^f/i(r) / dr'hnir')\U,{T - t')\' 
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/V^ ( 1 — 
cubic ' 



i d 
Ui dr 





(6) 



where ^ = and r = t/Ti in are dimensionless propagation axis and temporal de- 

lay (in the moving coordinates frame), respectively, Lqi is the dispersion length of 
the FW and Ti in is input pulse duration. Ui is the dimensionless envelope amplitude, 
-^1 ~ X]m=2 ^ m! ^'^^^ 7^"''" ( J^)"^ indicates the dimensionless dispersion operator for the 



i.e. iV4, = Lo.i/in^ 



FW, in which kj^\ujj) is the m-order derivative of kj at uj and Uj = cj^Ti^in is the di- 
mensionless frequencies, h^lr) is the nonlocal cascading response function whose normalized 
spectral function is hc{0,) = AA;(a;i)/AA;noniocai(^)5 ^^noniocai is called nonlocal phase mis- 
match which will be discussed later. /iR(r) is the material Raman response function. The 
dimensionless soliton numbers, A^^casc and A^cubic, are linked to the nonlinear refractive indices, 

ncasc and A^^ubic = -^D,i-^in^'"'cubic, where Jin is the peak intensity of the 
launched pump pulse, ncubic is the nonlinear refractive index induced by material cubic non- 
linearity which consists of electronic Kerr response nKerr,ei = (1 ~ /R)'^cubic and vibrational 
Raman response. 

4. Conditions of Cascading Nonlinear Soliton Pulse Compression 

In the limit of large phase-mismatch, cascading nonlinearities are produced through the 
coupling between the FW and the harmonics. Normally, ricasc dominates the cascading non- 
linearity and, particularly in SHG process, it can be expressed in the form of nonlinear 
refractive index [l]: (in unit of m^/W) 

(2) _ ^^l^eff (7) 

""case _ „9„2„ A v / 



eoC^nln2Ak 



6 



where Ui is the angular frequency of the FW, d^g = x*-^V2, and n2 are refractive indices 
of the FW and the SH, respectively. Equation ([t]) means that an equivalent Kerr nonlin- 
earity is generated by the quadratic nonlinearity, and the negative sign indicates that this 
nonlinearity shows an effect of self-defocusing with a positive Ak and vice versa. We remind 
the corresponding electronic Kerr nonlinearity nKerr,ei in the material, if the FW photon 
energy is sufficiently far away from the bandgap in the material, this nonlinearity is pos- 
itive (self-focusing) [24] and therefore forms a competing nonlinearity. The self-defocusing 
cascading nonlinear soliton pulse compression is therefore only supported if the Kerr-like 
cascading nonlinearity can counterbalance and exceed the material Kerr nonlinearity and 
form an overall effect of self-defocusing in the crystal. 

Hence, we get the basic condition of our pulse compression, i.e. ncasc > |?T-Kerr,ei| , and from 

(2) 

this arises a critical boundary of the phase-mismatch, Akc, for which ricasc + ^Kerr,ei = 0. An 
example is shown in Fig. [l](a, b), using MgOLN cut for Type-0 interaction. Within a broad 
wavelength region (950~1920 nm), the self-defocusing ncasc is stronger than the self-focusing 
^Kerr,ei while the dispcrsiou remains normal. In this range an overall effective self-defocusing 
nonlinearity can be found. 

Besides the basic condition of an effective self-defocusing nonlinearity, another condition, 
determined greatly by the nonlocal phase mismatch Afcnoniocai(fi), is necessary if ultrafast 
interaction is desired as it has a significant infiuence to the quality of pulse compression. 
One can namely derive a certain threshold of the phase- mismatch, Ak^^, below which the 
response of the nonlocal cascading response is resonant (non-stationary), and thus unsuitable 
for few-cycle interaction, while above it is non-resonant and has ultra-broadband (stationary), 
i.e. ideal for few-cycle interaction. Mathematically, when including complete dispersion the 
nonlocal phase mismatch can be written as [16]: 



AA;„oniocai(f^) = k2 {u2 + Q) - 2k^ (c^i) - Uk^^^ M = - ^^12^^ + Ak (8) 



m=2 



where (jJ2 is the frequency of the SH. k^p is the group velocity and di2 is the GVM between 
the FW and the SH, i.e. di2 = kP{oji) — k2\2uji). Essentially this equation describes the 
dispersion of the harmonic as observed from pump reference frame and the term Vtk^ {uji) 
is physically included to change the moving coordinate system from the the laboratory 
frame to the FW co-moving frame. A non-resonant cascading response is obtained when 
Afcnoniocai(^) 7^ for any (physical) frequency. The threshold phase- mismatch value Ak^r 
marks the boundary where some particular frequency inside the FW spectrum starts to 
experience nonlocal phase-matching according to this condition. This threshold can be ana- 
lytically obtained as d\2/2kf^ if higher order of the SH dispersion is neglected 14 



15 



or like 
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here, be numerically found. In Fig. [T|b) we show that in Type-0 MgOLN, the intrinsic Ak 
of the material is well within the stationary region. Therefore a clean pulse compression is 
expected, which can be explained as follows. 

Physically, in the non-stationary region, a nonlocal phase-matching condition will be ful- 
filled between a frequency in the side band of the generated harmonic and the FW center 
frequency. Due to the narrow bandwidth and large gain of this resonance, a strong oscillation 
occurs in the temporal profile of the harmonic pulse, resulting in degradation on the com- 
pression quality. On the other hand, if Ak stays above threshold in the stationary region, 
such a phase-matching condition will not be fulfilled and therefore the harmonic is merely 
slaved to the FW in the usual manner through A2 oc hcF[Uf] (as a part of the cascading 
channel necessary for generating a FW nonlinear phase shift). In this case the bandwidth of 
the cascading can be extremely high, even octave-spanning if the phase-mismatch is not too 
close to the resonance threshold Q , and in this condition few-cycle soliton formation with a 
clean temporal profile of the pulse is possible. 




Ai^^ (mm"') k (nm) 



Fig. 1. Compression window of MgOLN cut for Type-0 interaction, (a) Cas- 
cading quadratic nonlinearity (using ^33 = 25pm/V at 1064 nm and Miller 
scaling to other wavelengths) and the native electronic Kerr nonlinearity from 
the so-called two-band model [4||24]. (b) Critical boundary of producing an 
overall self-defocusing nonlinearity and the stationary/nonstationary thresh- 
old of the cascading nonlocal response function, (c) Radiation in the nonsta- 
tionary region, showing the resonant wavelength as a function of the effective 
phase-mismatch, with the pump located at 1560 nm. (d) Degradation on the 
compression window if using a 1*^* order QPM on MgOLN. 
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When short and broadband pulses interact through (shghtly) phase-mismatched harmonic 
generation, it is common consensus to interpret any resonances in the sidebands of the har- 
monic spectrum as simply being phase-matching points to a sideband of the FW spectrum, 
whose frequency detuning is such that the basic harmonic frequency locking is kept, in SHG 
i.e. ^2(1^2) 7^ 2^1 (wi), but k2{uj2 + VL) = 2ki{(jji + ^). However, this phase-matched sideband 
theory is clearly in contradiction with the nonlocal phase-matching theory here (Eq.([8])), 
which dictates that a resonance might occur if a sideband of the SH spectrum can become 
phase-matched with the FW spectrum taken at its center frequency. We now show that the 
nonlocal interpretation is the correct one by investigating the spectral behavior of harmonic 
generation in the non-stationary regime. To do that the phase-mismatch is reduced to below 
the resonant threshold, which can be achieved in Type-I interaction scheme by employing 
birefringent angle tuning close to the phase-matching point, or as we do here by employing 
quasi-phase-matching (QPM) technology in the Type-0 interaction scheme. Several simula- 
tions are done with different phase-mismatches and the resonant wavelengths are detected. 
Figure [T|c) shows the good agreement between the nonlocal theoretical prediction and the 
simulation results. Instead the phase-matched sideband theory does not give the correct 
description except very far from the resonant threshold predicted by the nonlocal theory. 

5. QPM-induced Influence and Self-steepening Performances 

QPM technology is well known as a tool to adjust phase- mismatch or even achieve an effective 
phase-matching. However, the price of using QPM technology is a reduced 2^^^^ order nonlinear 
constant due to the pre-factor in the modulation function. For example, is a factor ^ 
smaller if the 1^* order QPM is employed in the material. A graphical representation of this 
is that A/cc is reduced as shown in Fig. [l|^d), while the benefit is that the phase- mismatch 
can be flexibly tailored rather than being fixed (in Type-0 interaction scheme). We should 
therefore rewrite the cascading quadratic nonlinearity in a material with 1*^* order QPM 
structure as: 



(2) _ 2a;i(!4ff) . . 



where Afcgff = AA; ± ^ and A is the poling period. 

Since the reduction of d^^ is inevitable in the QPM structure, one needs to reduce A/cgfr to 
a level low enough so as to benefit from an increased cascading strength. However, doing this 
will have some other consequences. 1) The reduced A/cgfi will run the risk of getting into the 
resonant region. As is shown in Fig. [l|d), a break-even phase- mismatch value Afcbaiance lies 
where nf^^^ Qpjyj = n^Xc and it is mostly in the non-stationary region. Therefore in the par- 
ticular case studied here, Type-0 LN, QPM comes with the price of having a non-stationary 
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cascading nonlinearity except for longer wavelengths close to the zero-dispersion point. 2) 
The GVM- induced self-steepening effect of cascading will be enhanced as it is scaled by ^k~g. 
It will add to the system intrinsic self-steepening term and therefore give rise to stronger 
shock front on FW pulses. Moreover, with normal dispersion, self-steepening induced slow 
pulses will have chance to balance the Raman induced fast pulses. 

To analytically prove the latter consequence, we make use of the NLS-like equation. Besides 
the system's intrinsic self-steepening factor, the self-steepening induced by the cascading 
is revealed through the 1*** order response. In the weakly nonlocal limit, i.e. t/^(r — r') 
is assumed slowly varying compared with the response function. Hence, it can be Taylor 
expanded as Uf{T — r') ^ UI{t) — T'-^Uf^r). The nonlocal cascading response (the second 
term) in Eq.dol) therefore has: 



i d 
Cji dr 



f/*(r) 1 
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Coidr J ^ \ C02 dr 
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where = h'^{0) = (ii2/AA;((X'i)Ti in is defined as the response time of the cascading 1 
perturbation and iTc = J dr ■ t ■ /ic(t). It is obvious that a GVM-induced self-steepening 
term is produced which is scaled by Tc and will directly change the amplitude of the field 



envelope, resulting in shock front on pulses 11 



On the other hand, the case of the material Raman effects is quite different. Physically, 
material Raman effects are well understood to working only on the phase of the field envelope, 
chirping pulses, transferring energy from the high-frequency components to low-frequency 
components and resulting in red shifting in the spectrum. Mathematically, if we analogously 



expand the field intensity |f/i(r — r')| to the 1*^* order, the third term in Eq.(pj) has: 



"^cubic 



i d 
Ui dr 



cubic 



(1 - Jr) \U^{r)\'U,{r) + fnU^{r) / dr'h^{T')\U,{r - r')\^ 
i d 
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Ui dr 



\U^fU,~T^U,-^\U^\' 
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I 



+ HOT (11) 



where = f-^ j dr ■ t ■ hulr) is Raman response time. Note that the definitions of Tc and 
tr are different as the response time should be physically real-valued. Then, a real part is 
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produced as the l'^* order of the material Raman effects which directly works on the phase 
of the field envelope. 

Furthermore, the field envelope can be separated into a real-valued amplitude multiplying 
a phase term 12 , i.e. Ui = A{^,T)e'^'^^^''^\ In the dispersionless approximation, where the 



, i.e. Ui 

propagation length is much smaller than the dispersion length, the change on the amplitude 
and phase can be approximately obtained as: 



-Q^ = ( -. + 2^ciVcasc 
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(13) 



— --en— - ■ rv- ■ cuDic" 

where the effective soliton number, Nes = V^casc ~ ^cubic defined to quantify the overall 
self-defocusing nonlinearity. It is obvious that the cascading response time gives contributions 
to both the amplitude and phase change while the Raman response time only influences the 
phase change. 



Wilh disperaion: slow pulses (picfiicrar <0| 

Using QPM.it^=129Jmin"' 

No QPM, Ali=3l9.Smii 




(a) 




Fig. 2. Simulations of NLS-like equation on the balance of the fast/slow pulses. 
Dispersion properties are chosen from a LN crystal. The cascading nonlinearity 
is Ncasc = 2.60 and the material electronic Kerr nonlinearity is A^Kerr,ei = 1-96, 
the input pulse has a FWHM = 50 fs and a peak intensity /peak = lOOGW/cm^, 
located at 1560 nm. (a) shock front on pulses and slow pulses driven by dis- 
persion, (b) slow pulses in QPM crystal without Raman effects, in which 
Afceff = 129. 5mm^^. (c) fast pulses driven by 50% Raman effects in a bulk 
crystal with AA; = 319. 5mm'^. (d) balanced pulses generated in a QPM crys- 
tal with 50% Raman effects. The inset figures correspond to the spectrum evo- 
lutions, the temporal profiles and the spectra at the first compression point of 
the pulses. 



The self-steepening induced shock front on pulses is therefore understood to happen when 
the pre- factor on the right part of the amplitude equation (Eq.(12)) is not zero, see Fig. |2](a). 
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in which and N^asc play a role as tuning parameters. A similar conclusion was reported by 



Moses and Wise 11 but they did not get further into the amplitude and phase dynamics. 
Moreover, the shock front also causes asymmetry in pulse spectrum and, in the existence of 
the material dispersion, therefore generates fast/slow pulses. Fig. |2](a) also shows that for 
lower (effective) phase mismatch, slower pulses are generated. 



For the phase equation (Eq.(13)), we note that if a stable soliton is formed during the 



propagation, the phase chirp becomes negligible (the third term can be eliminated as ^ = 
for a stable soliton) and the Raman term turns to dominate (the second term), which contin- 
uously gives rise to red shifting in pulse spectrum. Then, with normal dispersion, fast pulses 
are generated due to the faster GV at the red-shifted wavelengths. 

Hence, we could expect a balance of the fast/slow pulses with the competing interaction of 
the cascading effects and the material Raman effects. Convincing numerical results are shown 
in Fig. |2](b-d) through direct simulations of Eq.([6]), in which the cascading nonlinearity and 
the material electronic Kerr nonlinearity are always kept identical. Without Raman effects, 
the cascading response time will be tuned by phase mismatch parameter and generates 
slow pulses, see Fig. |2](b); When the Raman term is switched on as 50% of the total cubic 
nonlinearity, fast pulses are driven as a counterbalance to the cascading effects. In a bulk 
crystal where the cascading response time is quite short due to the large Ak, eventually 
fast pulses are generated, see Fig. |2{c). When using a QPM crystal, however, the cascading 
response is enhanced and may form a balance with the Raman effects, which keeps the 
pulses around the zero delay position, see Fig. ^d). The oscillations split out of the pulses 
are dispersive waves (D-wave, also called Cherenkov wave), which are the result of a phase- 
matching condition of the linear part of the spectrum (the anomalous dispersion regime at 
longer wavelengths) to the soliton residing in the normal dispersion regime. 

6. Simulations and Discussions 

In this part, realistic simulation examples based on our master equations are shown. Few- 
cycle soliton pulses are generated within a short crystal length from multi-cycle pump pulses. 
Such a compressor has therefore great potential for improving the performance of femtosec- 
ond laser systems. 

The samples we take into consideration are BBO, MgOLN and PPMgOLN, which are all 
uniaxial crystals within the point group 3m. If these crystals are operated far from the band 
gap resonance, they are assumed to be under the Kleinman symmetry. Thus there are only 
a few non-zero elements in the nonlinear tensors, e.g. only 3 independent elements in 2"*^ 
order nonlinear tensor and 4 in the 3'''^ order nonlinear tensor. The nonlinear parameters 
used in the simulations here are listed in table [Tj The cubic nonlinear parameters of BBO 
were taken from a recent literature study [26], while those for MgOLN were taken from |;4|. 
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Note that the large Raman fraction of LN was found there by comparing simulations with 
experimental data, as no reliable absolute Raman nonlinear coefficients could be found in 
the literature. 

Table 1. Nonlinear parameters used for simulations 



Material 



BBO (1.03/im) 



MgOLN (1.3/im) 



PPMgOLN (1.3/im) 



2nd Q2~(-[gr nonlinear 
tensor elements 
(pm/V) 



d22 



2.20 15 
-0.04 
-0.04 



/ 
/ 

-23.50 



/ 
/ 

-23.50 



grd Qj-(^Qj- nonlinear 
tensor elements 
(pmVV^) 



Cll 

C33 
C16 

ClO 

/r 



550.00 
-1400.00 
120.00 
-22.00 




26 



26 



26 
26 



/ 

7300 

/ 
/ 

50% 



/ 

7300 

/ 
/ 

50% 



Interaction type: 



Type-I, Type-II 



Type-0 







/ ^ — 






































Type-0, QPM 



.LLTLTL 



6. A. Type-I BBO at 1030 nm: A Clean Few-Cycle Soliton Pulse Compressor 

A BBO crystal cut for Type-I interaction has a compression window for longer near-IR 
wavelengths (where the two compression conditions are both satisfied) which can be 
accessed by tuning the crystal rotation angle 9. Within this window, clean pulse compression 
is expected. In fact the successful generation of few-cycle compressed pulses at 1260 nm by 
Moses and Wise was done at such conditions. In Fig. |3] we investigate a similar scenario at 
1030nm where high energy sub-picosecond pulses can be generated by an ytterbium-doped 
laser system. While the phonon response of BBO has been measured 25 , its strength is 
generally believed to be very weak, so in the simulations we take /r = 0. Moreover, it is 
relevant to mention that BBO has normal FW dispersion below 1450 nm. 

The electric field of the compressed pulse is shown in Fig.[3]^a). After few centimeters prop- 
agation, clean and sub-3-cycle pulses are generated. Besides, the cascading response causes 
shock front on the pulses and generate slow pulses through normal dispersion, see Fig. |3|^b) 
showing the temporal evolution. In the spectral evolution Fig. [3]^c, d), the FW and several 
harmonics are marked and cascading nonlinearities are produced among them. A weak SH is 
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e-pol:IU l"(dB) 




-200 -100 100 200 1 2 3 4 5 6 7 



Fig. 3. Simulation of cascading nonlinear soliton pulse compression in 
BBO cut for Type-I interaction at 1030 nm. The rotation {9, (p) = 
(20.5°, 90.0°), under which Ak^ = rS.lmm'^ and Ah^ = 46.1mm-^ Na- 
tive phase-mismatch Ak = 55.4mm~^. The launched pulse has FWHM 
= 200 fs, peak intensity /i„ = lOOGW/cm^, niaL = -7.86 x lO'^VVW, 
^^Kerr,ei|o:ooo = 5.67 X lO'^VVW, A^'eff = 6.00. (a) The clcctrlc field of the 
stage soliton; (b) The temporal evolution of the FW; (c, d) Spectral evolutions 
of both o-pol and e-pol pulses. 
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generated through the phase-mismatched SHG process (Type-I, o+o:e) while a weaker third 
harmonic (TH) is generated through the phase-mismatched sum frequency generation (SFG) 
process (Type-II, o-|-e:o) as well as the phase-mismatched third harmonic generation (THG) 
process (0+0+0:0). Since the SHG process has the smallest phase-mismatch (~ 55mm~^) 
and moderate conversion efficiency compared with other processes, the cascading quadratic 
nonlinearity dominates the total self-defocusing nonlinearities. 

It is to our knowledge the ffist time that the simulation of the NWEFs including all the 
nonlinear processes (x'-^-' and x^^'') is presented and the electric field is directly modeled, 
which imply a generalized conception and simple expression. Moreover, we performed a 
comparison between the results of the NWEFs and the usual CWEs and found only a slight 
difference in the pulse's pedestal; this difference is caused by the THG that is included 
in the NWEFs, but since the TH is weak here only a slight difference is observed. Thus, 
the type I simulations in BBO traditionally done with the usual coupled-wave equations 
based on the SEWA |2||3||ll||T2|[T4}jT6) turns out to give a fairly accurate description of 
the dynamics in cascaded soliton compression. Nonetheless, the advantage of the NWEFs is 
that it automatically takes into account any type of multistep wave mixing and harmonic 
conversion, which cannot a priori be excluded. 

6.B. Type-0 MgOLN: Few- Cycle, Fast Pulses Dominated by the Material Raman Effects 

MgOLN cut for Type-0 interaction is also a good candidate for cascading nonlinear pulse 
compression due to its large ^33 element. The compression window has already been shown 
in Fig. [1} However, in MgOLN, it seems that the Raman effect is quite dominant |^, giving 
a large relative fraction /r, which causes fast pulses in time domain and red-shifting in 
spectrum, see Fig. |4j 

By launching a short 50 fs pump, 3-cycle pulses are generated within a short propagation 
length of around 4mm (Fig. |4|^a, b)) because MgOLN has a large GVD. It is highlighted 
that fast pulses are dramatically driven after the soliton formation in which the Raman 



chirping term starts to dominate the phase dynamics (Eq.(13)). Correspondingly, obvious 
red shifting in spectral evolution is observed and the D-wave becomes more blue shifted, see 
Fig. |4|^d), a consequence of the phase-matching law between them. Compared with Zhou's 
experiment [i], the pulse splitting is not prominent here as we use a lower soliton order and 
therefore reduce the effect of Raman-induced soliton fission. 

6.C. Type-0 PPMgOLN: Few-cycle, Zero Delay Pulses Generation with Enhanced D-wave 

To make a comparison with the bulk MgOLN crystal, pulse compression in a PPMgOLN 
with an effective phase-mismatch in the non-stationary region is shown in Fig. [5j Using the 
same pump, the compressed soliton has a single-cycle spike, but overall the soliton is not as 
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Fig. 4. Simulation of cascading nonlinear soliton pulse compression in 
MgOLN cut for Type-0 interaction at 1300 nm. The rotation [6, if) = 
(90°, 90°), i.e. X-cut, under which Ak^ = 890.4mm-i and Aksr = 283.3mm-^ 
Native phase-mismatch Ak = 501.5mm^^. The launched pulse has FWHM 
= 50 fs, peak intensity = lOOGW/cm^, n^iL = -40.05 x lO-^OmVW, 
nKerr,ei = 22.56 X lO-^Om^/W, ^A^isc - ^Kerr,ei = 2-05. (a) The electric field 
of the 1^* stage soliton; (b) The temporal evolution of the FW; (c) The spec- 
trum of the 1*^* stage soliton; (d) The spectral evolution of the e-pol pulse with 
obvious red shifting of the FW and blue shifting of the D-wave. 
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Fig. 5. Simulation of cascading nonlinear soliton pulse compression 
in PPMgOLN cut for Type-0 interaction at 1300 nm. The pol- 
ing period is 26 fim. The rotation {6, ip) = (90°, 90°), un- 
der which AfccQPM = 360.9mm~^ and Aksr = 283. 3mm^^. The effective 
phase-mismatch Akes = 259.8mm~^. The launched pulse has FWHM = 
50 fs, peak intensity 1;^ = lOOGW/cm^, n2L,QPM = -31-32 x lO-^OmVW, 

nxerr.ei = 22.56 X lO'^^m^/W, ^iVl,,QpM - iV2^„,ei = 1-45. (a) The electric 
field of the pulse; (b) The temporal evolution of the FW; (c) The spectrum of 
the pulse; (d) The spectral evolution of the pulse with resonant SH radiation 
and strong dispersive wave. 
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clean as Fig.|4]and has a stronger shock front. With a great reduction in the phase-mismatch, 
from 501.5mm~^ to 259.8mm~^ through a 26 fim QPM pitch, the 1*^* order response of the 
cascading (scaled by Tc) is enhanced which pulls the pulses back to the zero delay position, 
see Fig. |5]^b) showing the temporal evolution. In the spectral evolution Fig. ^d), Raman 
induced red shifting is well suppressed and the D-wave is evoked and enhanced around 
4.5/im. Besides, the self-steepening pre-factor is also enlarged, resulting the strong shock 
front on the pulses, see Fig. [sj^a). Meanwhile, in the non-stationary region, nonlocal phase 
matching is fulfilled so that a strong resonant peak emerges in the spectrum, see Fig. |5](c,d), 
which generates strong temporal oscillations and degrades the compression quality. 

7. Conclusions 

As a conclusion, in this paper we presented the generalized nonlinear wave equations in 
frequency domain (NWEFs), which for the first time include polarization coupling, second 
and third order anisotropic nonlinear it ies and delayed Raman nonlinear effects. Any of the 
possible phase-matching conditions (Type-0, -I and -II) are easily modeled, and by including 
the anisotropy of the nonlinearities, the different wave-mixing possibilities (3- and 4-wave 
mixing, including third-harmonic generation and parametric up- and down-conversion) are 
automatically included properly. We then used the NWEFs as a platform for investigating 
ultrafast cascaded second-harmonic generation, in which soliton compression to few-cycle 
duration is possible in short nonlinear frequency conversion crystals. 

We first reduced the NWEFs to a single NLS-like equation, which highlights self-steepening 
effects as well as the competition between cascaded nonlocal effects and cubic SPM and 
Raman effects. Using this reduced equation we discussed the conditions for observing optimal 
pulse compression. It is common to interpret the resonant SH peak, which can be observed 
in cascading when the phase mismatch is small, as phase-matching of the spectral sidebands. 
We showed that this is wrong and that it can instead be accurately described by the nonlocal 
theory as a phase-matching between an offset SH frequency with the FW center frequency. 

The reduced NLS-like equation was also used to investigate the competition between self- 
steepening, Kerr self-focusing SPM, and Raman effects on one side, and on the other side 
cascaded contributions. Besides the self-defocusing SPM term, these cascading terms to first 
order contribute both to self-steepening-like terms (creating a pulse front shock) and also to 
higher-order phase contributions similar to the Raman effect. Since the cascading terms are 
tunable through the sizes and signs of the phase-mismatch and the GVM parameters, we 
showed that the slow/fast pulses usually associated with self-steepening and Raman effects 
could be balanced out by cascading: The Raman red-shift of the pulse gave fast solitons due to 
the faster GV at red-shifted wavelengths (normal dispersion), while GVM-induced cascading 
tends to slow pulses down. Therefore, pulses where the two delay effects are balanced are 
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expected and we showed simulations where these effects appeared. In fact, we even found 
cases where the strong spectral red-shift from Raman was completely suppressed by strong 
cascading contributions. 

Finally, we showed a series of full simulations of the NWEFs, where few-cycle, clean 
soliton pulse compression was demonstrated. In BBO the Type-I simulations showed spectral 
components related to third-harmonic generation and sum-frequency mixing, but we found 
that the results agreed well with previous simulations of the coupled SEWA equations. In 
MgOLN a strong Raman contribution was investigated in the Type-0 configuration, and we 
showed that in a bulk crystal a few-cycle soliton forms after a few millimeters, but also that 
due to the dominating Raman effect the soliton is accelerated. We also performed the same 
simulation with the same pump but in a QPM (through periodic poling) MgOLN, where the 
phase mismatch was greatly reduced. In this case the cascading response is no longer optimal 
as it becomes spectrally resonant. The QPM case gave more characteristic phenomena in the 
spectrum but the temporal profile was not so clean. In addition, the soliton was no longer 
accelerated by Raman effects, and spectrally the Raman red-shifting was cancelled. This is 
a consequence of the increased cascading self-steepening from QPM. 
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